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ABSTRACT 

It is well known that newly formed planetary systems undergo processes of orbital reconfiguration and planetary migration. As a 
result, planets or protoplanetary objects may accrete onto the central star, being fused and mixed into its external layers. If the 
accreted mass is sufficiently high and the star has a sufficiently thin convective envelope, such events may result in a modification 
of the chemical composition of the stellar photosphere in an observable way, enhancing it with elements that were abundant in the 
accreted mass. The recent Gaia-ESO Survey observations of the 10-20 Myr old Gamma Velorum cluster have enabled identifying a 
star that is significantly enriched in iron with respect to other cluster members. In this Letter we further investigate the abundance 
pattern of this star, showing that its abundance anomaly is not limited to iron, but is also present in the refractory elements, whose 
overabundances are correlated with the condensation temperature. This finding strongly supports the hypothesis of a recent accretion 
of rocky material. 

Keywords. Stars: abundances - Stars: chemically peculiar-Stars: pre-main sequence-open clusters and associations: individual: 
Gamma Velorum - Techniques: spectroscopic - Planet-star interactions 


1. Introduction 


Open clusters are groups of stars formed from the same nebula. 
For this reason, members of the same cluster are expected to 
share distance, age, kinematics, and chemical content. The latter 
reflects the initial composition of the cloud that gave birth to the 
cluster itself. 

Recently, a few cases of chemically anomalous stars have 
been observed in open clusters; these can be explained through 
episodes of planet engulfment or the significant accretion of 
rocky material (Ashwell et al. 2005; Spina et al. 2014). Stars 
are subject to occasional accretion events during their lifetime 
through different processes. In addition to the well-known phase 
of gaseous accretion that characterizes newly born stars during 
their first Myr of evolution, infall of planets or planetesimals 
onto the central star can represent another important mechanism 
that is driven by orbital decay following interactions with other 
massive bodies (e.g., Weidenschilling & Marzari 1996; Kley & 
Nelson 2012). The most favorable epoch for such events extends 
approximately from the end of the main protostellar accretion 
phase to the completion of pre-main-sequence (PMS) contrac¬ 
tion, when newly formed planets are still clearing their orbits of 
dust and rocky bodies. A major consequence of the ingestion of 
(rocky) planetary material is the possible enhancement of pho- 
tospheric metallicity: after penetrating the stellar atmosphere, 
rocky masses are rapidly dissolved and mixed with the ambient 
matter. If the accreting star has a thin convective zone (CZ), the 
polluting planetary material is not too diluted within the predom¬ 


inantly hydrogen gas and can produce a significant increase of 
the atmospheric metallicity (Laughlin & Adams 1997; Gonzalez 
2006; Schuler et al. 2011; Mack et al. 2014). Such an enhance¬ 
ment might be detected through high-resolution spectroscopic 
observations. 

The Gamma Vel cluster is a young (10-20 Myr; Jeffries 
et al. 2014) open cluster that has recently been observed by the 
Gaia-ESO Survey (Gilmore et al. 2012; Randich et al. 2013). 
As we will show, the age of this cluster offers the rare op¬ 
portunity to investigate the impact that accretion events may 
have on the composition of the stellar atmosphere. Recently, 
Spina et al. (2014) determined the cluster mean metallicity 
<[Ee/H]>= -0.057+0.018 dex and found that one of the mem¬ 
bers, 2MASS J08095427-4721419 (hereafter, #52 as in Spina 
et al. 2014) has an iron abundance 2cr greater than that of the 
cluster. The star also has an IR-excess at 24 pm (Hernandez 
et al. 2008) that is due to a long-lived debris disk that may have 
already formed planetary-mass objects. Assuming that #52 is 
genuinely metal-richer than other cluster members. Spina et al. 
(2014) argued that the overabundance might be produced by the 
ingestion of ~60 Me of rocky material. Here, we further investi¬ 
gate the peculiarity of #52 by comparing its chemical pattern to 
that of a similar star of the same cluster. We show that the ac¬ 
cretion scenario is confirmed by the overabundance of refractory 
elements found in the atmosphere of #52. 


2. Dataset 


* Based on observations made with the ESOWLT, at Paranal Ob¬ 
servatory, under program 188.B-3002 (The Gaia-ESO Public Spectro¬ 
scopic Survey) 


Gaia-ESO observations of Gamma Vel have been described by 
Spina et al. (2014) and Jeffries et al. (2014). In this paper we use 
UVES data, in particular, the results from the second analysis 
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cycle that are part of the second internal data release of the Gaia- 
ESO Survey (GES) in July 2014 (GESviDR2). Our analysis is 
based on the recommended abundances delivered by two work¬ 
ing groups (WGs); one in charge of the UVES spectra of EGK 
type stars (i.e., WGl 1), and one dealing with the spectra of PMS 
stars (i.e., WG12). Details of the analyses performed by these 
WGs are reported by Smiljanic et al. (2014) and Lanzafame et al. 
(2015). In both working groups the spectrum analysis is carried 
out by different nodes, whose results are then combined to pro¬ 
vide a set of recommended parameters. In addition to the rec¬ 
ommended abundances and stellar parameters provided by the 
GES, in this study we also exploit the abundances produced by 
different nodes of the consortium so as to explore the reliability 
of the analysis, as shown below. 


3. Chemical anomaly of #52 

Abundances of a variety of elements (Na, Mg, Al, Si, Ca, Sc, 
Ti, V, Cr, Mn, Ee, Co, Ni, and Cu) for two members of Gamma 
Vel, the metal-rich star #52 and 2MASS J08093304-4737066 
(hereafter, #45 as in Spina et al. 2014), have been determined 
in GESviDR2 by four nodes. Of all the cluster members, only 
these two stars are sufficiently slow rotators for a detailed chem¬ 
ical analysis, thus individual abundances could not be derived 
for other stars in Gamma Vel. However, we stress that star #45 
has a metallicity similar to that of the other cluster members, 
hence it is probably representative of the cluster mean composi¬ 
tion. We compare the chemical abundances derived for the two 
stars to gain more insight into the chemical pattern of #52. The 
spectra have a similar S/N (127 and 140, respectively), and from 
them, the following atmospheric parameters have been derived 
and released in GESviDR2: star #45 has an effective tempera¬ 
ture of 5614+80 K, a surface gravity of 4.22+0.12 dex, and a 
microturbulence of 1.85+0.08 km/s, while #52 has 5864+57 K, 
4.41+0.11 dex, and 1.89+0.10 km/s. 

Eor each element X, we considered the nodes that pro¬ 
duced a A(X)=logVx/logVH+12 value for both stars. Then, 
for each element, we computed the abundance differences be¬ 
tween #52 and #45 determined independently by each node: 
AA(X)=A(X)# 52 -A(X)# 45 ; these differential abundances are free 
from the systematics or zero-points that can affect the analysis 
of each node. The elemental abundance values of #52 recom¬ 
mended by WG12 (A(X)) and the AA(X) derived by each node 
are listed in Table 1. The differential abundances provided by the 
nodes show a satisfactory agreement for many elements, such as 
Si, Sc, and <Ti>, whose standard deviations around the mean 
are 0.05, 0.03, and 0.06 dex, respectively. However, for other 
elements such as Na, Ca, or Ni, one measurement appears to 
strikingly disagree with the others of the same element. To pro¬ 
vide a single AA(X) for each element, we therefore opted for the 
median of all the values since it is a robust measure of the central 
tendency of the differential abundances given for each element. 
Thus, in the last column of Table 1, we list the median values AA 
of all the elements together with the median absolute deviation 
(MAD). The uncertainties are typically <0.10 dex, with only one 
exception: that for copper. Its large error (i.e., 0.41 dex) is due 
to highly discordant determinations provided by two nodes. 

Using the IRAE splot task, we carefully measured the equiv¬ 
alent widths (EWs) of the Cu lines at 5105.52 A in the two spec¬ 
tra, obtaining 101 mA for #52 and 113 for #45. The expected 
uncertainties in the measured EWs can be estimated from the 
Cayrel (1988) formula, which predicts an accuracy of 1mA for 
these lines. We note that this formula neglects the uncertainty 


Table 1. Differences based on Gaia-ESO node analyses. 


Element 

A(X) 

AAi 

AA2 

AAs 

AA4 

AA+MAD 

Nal 

6.24 

0.62 

-0.06 

-0.16 


-0.06+0.10 

Mgl 

7.74 


0.33 


0.12 

0.23+0.11 

All 

6.66 


0.19 


0.06 

0.13+0.07 

Sil 

7.57±0.14 

0.12 

0.08 

0.04 

0.18 

0.10+0.04 

Cal 

6.48±0.14 

-0.26 

0.20 


0.19 

0.19+0.01 

Sell 

3.4!±0.20 

0.24 



0.29 

0.27+0.03 

Til 

5.16±0.27 

0.05 

-0.01 

0.13 

0.07 


Till 

5.31±0.36 

0.07 

0.39 

0.07 

0.30 


<Ti> 


0.06 

0.19 

0.10 

0.19 

0.15+0.05 

VI 

4.22±0.23 

0.27 


0.15 

0.01 

0.15+0.12 

CrI 

5.8±0.10 

0.13 

-0.02 

0.18 

0.08 

0.09+0.05 

Mnl 

5.84±0.21 

0.21 

0.06 

0.05 

0.15 

0.11+0.05 

Fel 

7.55±0.12 

0.08 

0.37 

0.12 

0.22 

0.17+0.07 

Col 

5.16±0.33 

0.07 

0.15 

0.11 

-0.09 

0.09+0.04 

Nil 

6.28±0.14 

0.15 

0.12 

0.42 

0.11 

0.14+0.02 

Cul 


-0.40 



0.41 

0.01+0.41 


a) <Ti>=(log Nn i + log A'r, ;/)/2 

b) The Gaia-ESO nodes are 1-EPINARBO, 2-LUMBA, 3-Vilnius, and 4-Aicetri. 


Table 2. Differences based on a direct analysis of the Gaia-ESO spectra 


Feature 

EW#52 

EW#45 

AA 

Cu (5105.52 A) 

101+1 mA 

113+1 mA 

0.00+0.04 

Zn (4810.56 A) 

92+1 mA 

93+1 mA 

+0.02+0.05 


due to the continuum setting, but this contribution is minimal 
in a differential analysis if adopting a similar assumption on the 
continuum location in both the spectra. Based on these EWs, 
assuming the stellar parameters mentioned above and the anal¬ 
ysis tools used by the GES for the DR2 analysis (i.e., atomic 
parameters by Heiter et al. 2014 and MARCS stellar models; 
Gustafsson et al. 2008), we obtained a differential abundance 
AA(Cu)=0.00+0.04 dex. This error takes into account the un¬ 
certainties in the atmospheric parameters of the two stars and 
has been evaluated as follows: we determined the abundance er¬ 
rors for #52 and #45 as in Magrini et al. (2013), then we adopted 
as final uncertainty the maximum difference of the abundance 
values considering their error bars. We neglected the uncertainty 
due to the EW measurements. A similar check has been per¬ 
formed on the zinc lines at 4810.56 A: their EWs are 92 mA 
for #52 and 93 for #45. This resulted in a differential abundance 
AA(Zn)= +0.02+0.05 dex. As for Cu, the uncertainty on the 
EWs (i.e., ~1 mA) are negligible. These differential abundances 
are listed in Table 2. By repeating the procedure for four Sc II 
lines, we verified that our method provides differential abun¬ 
dances that are consistent with those listed in Table 1. 

The median values A A are plotted in Eig. 1 as a function 
of the elemental condensation temperature (Tcond)- Eor the lat¬ 
ter we have adopted the values of Lodders (2003). We also in¬ 
clude the values for Cu and Zn from Table 2 (triangles) that are 
based on the direct measurements of the lines observed in the 
two stars. It is possible to separate all the elements plotted in 
the graph into three illustrative classes: the volatiles would be 
the elements with Tc„,id<i^00 K (i.e., Zn, Na, and Cu), moder¬ 
ately refractory elements those with 1100<7’con£/<1500 K (i.e., 
Mn, Cr, Si, Mg, Ee, Co, and Ni), and the super refractory ele¬ 
ments those with 7’co„£/>1500 K (i.e., Ca, Ti, Al, and Sc). Eor 
each of these classes we plot in Eig. 1 the mean and standard de¬ 
viation of the differential abundances AA weighted on the corre¬ 
sponding MADs (dashed lines and oblique lines): 0.00+0.03 dex 
for the volatiles, 0.127+0.014 dex for the moderately refractory 
and 0.196+0.009 for the super refractory elements. Thus, even 
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Fig. 1. Differential abundances of #52 and #45 as a function of the 
condensation temperature. The plotted AA values are those reported in 
Table 1. The triangles are from Table 2. The condensation temperatures 
are taken from Lodders (2003). The horizontal dashed lines represent 
the weighted means of the differential abundances within each class of 
elements (i.e., volatiles, and moderately and super refractory elements). 
The bands represent the weighted standard deviation of the mean. 

though open clusters are commonly found to be chemically ho¬ 
mogeneous (e.g., De Silva et al. 2006), the plot shows that all 
the refractory elements (7’co„rf>l 100 K) are enhanced in #52 by 
>0.10 dex. On the other hand, this overabundance is not ob¬ 
served in the volatiles (Zn, Na, and Cu). In addition to this, the 
super refractory elements are more enhanced than the moder¬ 
ately refractories. Moreover, the most refractory element (Sc) 
is the most overabundant. Therefore we conclude that the over- 
metallicity of star #52 with respect to #45 is only limited to the 
elements with 7’co„£/>1100 K and that the overabundance of a 
given element is correlated with its condensation temperature. 

4. Enrichment by planet engulfment 

This study represents a significant improvement on the analysis 
of the anomalous Fe abundance of #52 performed by Spina et al. 
(2014) since it is based on a subsequent internal GES data re¬ 
lease with abundances of an extended set of elements obtained 
by four independent groups. This allows us to place the peculiar 
chemical abundances on a firmer basis and to investigate now the 
possible origin of this anomaly, taking advantage of the youth of 
Gamma Vel, whose age corresponds to the most active and dy¬ 
namical phase of planet formation and early evolution. 

Interestingly, the trend between the overabundance of a given 
element and its Tcond shown in Fig. 1 mirrors the composition 
seen both in planetary material (Chambers 2010) and in inter¬ 
stellar dust (e.g., Spitzer 1978). This suggests that the abun¬ 
dance pattern observed in #52 originated from the enrichment 
subsequent to the ingestion of planetary-mass-sized rocky ob¬ 
jects. Because of their higher condensation temperature, the re¬ 
fractories are thought to be the main components of the solids 
that accrete onto planets or planetesimals. 

Several studies have demonstrated that episodes of planetary 
material engulfment must be frequent in the highly dynamical 
environment of young systems. For example, high-precision ra¬ 
dial velocity surveys have revealed that about 20-30% of solar- 
type stars have low-mass planets (1-30 Me) with orbital periods 


shorter than 50 days (Mayor & Udry 2008; Howard et al. 2010). 
Several numerical simulations indicate that a large number of 
terrestrial-planet embryos can form within the innermost regions 
(<1 AU) of a planetary system on timescales shorter than 10 Myr 
(e.g.. Chambers 2001; Ida & Lin 2008; Schlichting 2014). Mu¬ 
tual interactions and tidal perturbations of the host star can eas¬ 
ily induce their orbital decay and subsequent accretion onto the 
central star (e.g., Zhou 2010; Ida &. Lin 2008; Raymond et al. 
2011). Similarly, the high frequency of Jupiter-type planets with 
extremely small orbits clearly points toward efficient planet mi¬ 
gration from their birth sites at large distances. During this pro¬ 
cess, a giant planet can induce other planets to move into un¬ 
stable orbits (Zhou & Lin 2008). Such catastrophic events are 
expected to occur more frequently during the early stages of the 
evolution of planetary systems (10-100 Myr), when the chaotic 
growth of planetesimals into planets is still active (e.g., Naga- 
sawa et al. 2007, and references therein). 

With an age of ~ 10-20 Myr, the members of Gamma Vel are 
prime candidates for studying differential enrichment due to ac¬ 
cretion of planetary material. These stars are old enough to host 
planets or planetesimals that interact with the remnants of the de¬ 
bris disk, but are not too evolved to lose the evidence of the pos¬ 
sible modification of their surface abundance. We know that any 
induced additional metallicity is bound to decrease as a result 
of several physical processes operating on timescales of several 
tens to hundreds of Myr. One of the most effective mechanism 
for the removal of chemical anomalies is the so-called thermoha¬ 
line convection (e.g., Vauclair 2004), in which the heavy material 
accumulated in the external layers of a star produces an unstable 
weight gradient that triggers the onset of convective instabilities 
that rapidly sink the overmetallic matter into the radiative stellar 
interior. Theoretical studies predict that this process would be 
able to halve the induced overmetallicity in ~50 Myr (Theado & 
Vauclair 2012). 

The second aspect to consider is that the mass enclosed in the 
CZ of a star is a critical parameter for determining the amount of 
chemical alteration induced by the accreted planetary object. It 
is well known that, during the early phase of contraction, solar- 
type and intermediate-mass stars undergo a process of internal 
readjustment in which the extended inner CZ retreats toward the 
surface (Palla & Stabler 1993). Using the PMS models of Siess 
et al. (2000), we illustrate in Fig. 2 the variation in CZ mass 
(solid line) during the first 30 Myr of the evolution of a 1.3 Mq 
star, corresponding to the mass estimated for #52 from its T'eff, 
V mag, and assuming the Siess tracks. In this case, the extent 
of the CZ has shrunk to about 5% of the total mass after about 
15 Myr and disappears (with the exception of a thin subpho- 
tospheric layer) in another few million years. The figure also 
shows the time variation of the photospheric iron content of the 
star (cf. dashed lines) resulting from the accretion of objects with 
different masses during the PMS contraction phase: 1, 6, and 20 
Me of pure iron. The latter value can be considered as an upper 
limit that corresponds to the total amount of iron contained in 
all the planets of our solar system, assuming the elemental abun¬ 
dance distribution of the meteoritic Cl chondrites reported by 
Lodders (2003). As we can see, a fixed accreted mass of rocks 
can produce extremely different iron enhancements depending 
on the stellar age owing to the rapidly evolving decrease of the 
CZ. 

The age of #52 is estimated to be ~15 Myr (vertical dotted 
line in Fig. 2; Spina et al. 2014), a value consistent with that 
of Gamma Vel. From the figure we note that at this evolution¬ 
ary stage the star has a sufficiently thin CZ that the accretion 
of 6 Me of pure iron would be enough to yield the observed 


Article number, page 3 of 5 








iron overabundance (shown by the horizontal dotted line). This 
amount corresponds to ~30 Me of rocky material having the 
chemical composition of Cl chondrites (Lodders 2003). On the 
other hand, if #52 were slightly younger (for instance, ~ 10 Myr), 
the ingestion of all the iron present in the solar system, would 
have produced a A[Fe/H] <0.08 dex, which is lower than ob¬ 
served. Finally, we see that the ingestion of 1 Me of pure iron 
can also cause the observed variation of [Fe/H] if the stellar age 
is greater than ~17 Myr. Owing to the uncertainty in the stellar 
age and mass and in the models of PMS evolution, we cannot be 
sure of the exact values of all the parameters, but the exercise is 
useful to appreciate that to match the observed overabundance, 
the accretion episode must have occurred within the last ~5 Myr. 

The PMS models also indicate that, at the current age of 
Gamma Vel, the solar-type members are in the process of sig¬ 
nificantly reducing the extent of their CZs to the point where 
they become thin enough to show some evidence of chemical 
pollution. Therefore, we may expect that in addition to #52, 
other stars may have undergone similar episodes of planet en- 
gulfment in the past that substantially enriched their chemical 
composition. Flowever, such hypothetical metal-rich stars could 
be rare objects because the accretion events may have occurred 
more frequently during the earliest stages (<10 Myr) of plan¬ 
etary evolution when the stars still had extended CZs that di¬ 
luted the chemical contamination. The steep dependence on age 
of the retreat of the convection zone is another factor that very 
likely limits the probability of discovering anomalous stars in 
this or other young clusters. It is also conceivable that other 
planetary engulfment events might be triggered within the first 
several hundreds of Myr, during the physical reconfiguration of 
the planetary system architecture. However, the additional effect 
of thermohaline convection in reducing or eliminating any sign 
of enrichment in less than 100 Myr suggests that the only suit¬ 
able clusters in which to find such chemically peculiar objects 
are those with ages of between ten and a few hundreds of Myr. 


5. Star #52 in the context of anomalously enriched 
stars 

Chemical patterns similar to that of #52 have been found in stars 
of widely different ages. However, different interpretations have 
been proposed to explain these chemical peculiarities. 

Melendez et al. (2009) and Ramfrez et al. (2009) have found 
that the Sun is depleted of refractory elements with respect to 
a sample of solar twins. This has been interpreted as evidence 
of the condensation of refractory elements into planets that oc¬ 
curred in the solar protoplanetary disk, but not around the ma¬ 
jority of solar twins. The material locked up into planets did not 
fall onto the Sun during the main accretion phase and therefore 
did not pollute the solar CZ. Conversely, the solar twins that did 
not form planets were subject to the accretion of the disk ma¬ 
terial that enhanced the abundance of refractories relative to the 
Sun. Tucci Maia et al. (2014) have provided further support to 
this picture with the observation of a lack of refractory elements 
in the planet-hosting star 16 Cyg B with respect to its compan¬ 
ion, 16 Cyg A. The abundance pattern found by Tucci-Maia et 
al. is remarkably similar to that shown in Fig. 1. This suggests 
that the observed chemical anomalies have a similar origin in the 
formation of rocky bodies. However, we note that the overabun¬ 
dance of refractories found in #52 is considerably larger than 
that of 16 Cyg A. As we have shown before, the sensitivity of 
the enhancement to the size of the CZ at the time of accretion 
can account for the observed difference. Considering the age of 



Age (Myr) 


Fig. 2. Mass contained in the CZ of a 1.3 Mq star as a function of 
its age during the PMS contraction phase (solid line). The right scale 
shows the iron abundance level resulting from accretion. The plot as¬ 
sumes the Siess et al. (2000) models for pre-main-sequence stars with 
solar metallicity. The rise in iron abundance caused by the ingestion of 
1, 6, and 20 Me of pure iron (dashed lines, from right to left) is highly 
dependent on the thickness of the CZ and, consequently, on the stellar 
age. The dotted vertical line is located at the age estimated for #52, and 
the dotted horizontal line refers to the overabundance of #52 relative to 
Gamma Vel. 

16 Cyg A (7.3 Gyr), other mixing processes may have had time 
to reduce the initial enhancement. 

The pollution from self-enrichment of pure rocky material 
has also been proposed to explain the overabundance of refrac¬ 
tory elements found by Ashwell et al. (2005) in the 600 Myr 
old J37 F-type star with respect to the other members of 
the NGC 6633 cluster. Notably, the observed enhancement, 
A[Fe/H]=0.85, is conspicuous and higher than that found in #52, 
which has a much thicker CZ. 

Finally, in a sample of 148 solar-type stars Adibekyan et al. 
(2014) have found a clear correlation between abundances and 
Tcond similar to that found in the present study. However, they 
suggested that this is more likely related to stellar age (e.g., 
old stars have smaller abundances of refractory elements re¬ 
spect with volatiles) or to Galactocentric distances. Our analysis, 
based on the comparison of stars of the same cluster, indicates 
that the peculiar overabundances of refractories found in #52 can 
be reasonably explained through pollution processes that are due 
to external causes, such as accretion of rocky material. 

6. Perspectives 

The recent discovery of several stars with a distinctive chemi¬ 
cal abundance pattern showing an overabundance of refractory 
elements over the volatile ones and the correlation with the con¬ 
densation temperature has opened the observational study of ac¬ 
cretion events of rocky bodies (planets, planetesimals) onto their 
parent stars. The study of the anomalous star #52 of Gamma Vel 
highlights the importance of studying young systems when the 
dynamical interaction of forming planets is at a peak and other 
physical mechanisms operating on longer timescales have not 
yet had time to wash out possible chemical features. More ac¬ 
curate abundance determinations of both refractory and volatile 
elements can constrain the nature and the composition of the ob¬ 
jects accreted onto the star. Finally, observations of G- and F- 
type stars in open clusters with ages of between ten and a few 
hundreds Myr will provide the opportunity to investigate the fre- 
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quency of these chemically polluted stars, the possible causes 
that trigger the events of planet engulfment (e.g., the presence of 
a debris disk or a stellar companion), and the fate of the rocky 
material accreted onto the star. 
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